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Abstract. 

The compressive strength of unidirectional fiber composites in the form of Kevlar yarn 
with a thin outer layer of graphene was investigated and modeled. Such fiber structure 
may be fabricated by using a strong chemical bond between Kevlar yarn and graphene 
sheets. Chemical functionalization of graphene and Kevlar may achieved by modification 
of appropriate surface-bound functional (e.g., carboxylic acid) groups on their surfaces. 
In this report we studied clastic response to unidirectional in-plane applied load with 
load peaks along the diameter. The 2D linear elasticity model predicts that significant 
strengthening occurs when graphene outer layer radius is about 4 % of kevlar yarn radius. 
The polymer chains of Kevlar are linked into locally planar structure by hydrogen bonds 
across the chains, with transversal strength considerably weaker than longitudinal one. 
This suggests that introducing outer enveloping layer of graphene, linked to polymer 
chains by strong chemical bonds may significantly strengthen Kevlar fiber with respect 
to transversal deformations. 



PACS numbers: 62.20.-x, 68.35.Gy, 81.40.-z 
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1. Introduction 

Reinforced polymer fibers have been studied intensively in search of superior fabrics for 
a broad range of applications including bullet-proof vests, protective clothing, and high- 
performance composites for aircraft and automotive industries [U, [2] . For most structural 
materials, the compressive strength is much greater than the tensile strength. Fiber- 
reinforced composites are among a few materials that exhibit greater tensile strength 
than compressive strength. This behavior results from fiber micro buckling compressive 
failure mechanism in fiber composites. Static compressive strengths of unidirectional 
composites have been studied by several research groups in the past few decades [3-6] . 

A major step forward in the protection of personnel against ballistic threats has 
been achieved by introduction of Kevlar synthetic fibers and textiles. The most common 
approaches for the fabrication of reinforced composite fibers have been melt processing 
and solution coagulation spinning [7J [8] Kevlar fibers are produced by wet-spinning from 
sulfuric acid solutions [9J. While in acid solution, Kevlar form liquid crystalline phase, 
with rod-like polymer chains well-aligned parallel to each other, and bonded by hydrogen 
bonds across to form strong planar sheets [10]. Such sheets are stacked radially to 
form Kevlar fiber. As a result, the fiber is much stronger with respect to longitudinal 
deformations than to transversal deformations. 

On the other hand, the carbon nanotubes (CN) and graphene sheets (GS) have 
been proposed as one of the most promising additives for the fabrication of ultra-strong 
polymer composites due to their advanced mechanical properties [TTJ [T2] . It is well-known 
that carbon nanotubes and graphene sheets have Young's modulus and tensile strength 
above 1 TPa and 60 GPa, [13-16] respectively, while their densities can be as low as 
1.3 g/cm 3 . They may exist as single layered or multi-layer structures. It is possible 
to harness the multifunctional properties of graphene sheets and design novel class of 
advanced composites with superior mechanical and electric performance p21 HE] Thus, 
various polymeric-CN and polymeric-GS composites and their chemical functionalisations 
have become a subject of intensive research and technological development over the last 
years [19j [20] . Recently, a significant mechanical enhancement of Kevlar fibers has been 
demonstrated by incorporation of surface modified carbon nanotubes achieved by using 
N-methylpyrrolidone (NMP) solvent [2T1 [22]. The maximum values observed were for 
Young's modulus 115-207 GPa; yield strength 4.7-5.9 GPa; strain at break 4.0-5.4 %; and 
toughness 63-99 J/g [21]. 

In this paper we present the 2D linear elasticity model to predict the mechanical 
properties of the Kevlar fibers reinforced by depositing an outer shell of graphene, rigidly 
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attached to Kevlar fiber by chemical bonding. The schematic diagram of the kevlar- 
graphene fibers is shown in Figure l,(a). Graphene outer layer can be formed by soaking 
kevlar single fiber with graphene dispersed in N-methylpyrrolidone (NMP) solvent. By 
using spinning technique it will be possible to create a composite structure by depositing 
an outer layer of graphene enveloping the Kevlar fiber. This lead us to hypothesize that 
graphene could be incorporated into swelled Kevlar fibers by soaking fibers in a dispersion 
of graphene in NMP. Chemical functionalization of graphene and Kevlar may be achieved 
by modification of appropriate surface-bound functional (e.g., carboxylic acid) groups 
on their surfaces. We anticipate that outer graphene layer enveloping Kevlar fiber may 
stabilize radially stacked hydrogen-bonded planar sheets forming Kevlar fiber with respect 
to transversal deformations. 

The main goal of this paper is to model the mechanical behavior of kevlar-graphene 
fibers and to predict the optimal thickness of the outer layer of graphene for strengthening. 
We show that within 2d linear elasticity model, and taking realistic Young modulus and 
Poisson ratios for transversal deformations, reinforcement starts to occur when the outer 
hard shell radius is about 4 % of kevlar fiber radius (corresponding to about 8 % weight). 
We assume that applied load has a narrow load peak along the diameter, which we may 
take to be the Y axes; typical displacements are shown in Figure l,(b). 

Our modeling assumes that graphene flakes are tied up together by covalent bonds, 
which may be achieved by graphene functionalization. In the presence of strong interlayer 
covalent bonds, we assume that bonded graphene flakes will have high Young modulus 
and tensile strength in all directions, and we model it by isotropic continuous model. Such 
model predicts theoretical limit on strengthening which may be achieved. 

Commercially available Kevlar fiber has radius of about 6 microns. In that case, 
we will need about 240 nanometers of functionalized multilayer graphene for substantial 
two-fold yield strength improvement. At such scale, continuous elasticity model may be 
used. 

2. Explicit solution 

We assume that the enveloping graphene layer and Kevlar fiber are rigidly attached by 
chemical bonds on their common boundary, and that the load is perpendicular to the 
fibers. 

In this section we derive an explicit solution of planar linear elasticity equations, 
modeling a circular fiber enveloped by a rigidly attached to it ring of different stronger 
material. 
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In section 3, we use our explicit solution to compute stress tensor, both in the outer 
graphene layer and inside Kevlar fiber, and to analyze increase of the yield strength of 
the composite system, as a function of reinforcing layer thickness. 

Our explicit solution is given in terms of the Airy stress function U. Airy stress 
function is a function from which the stress (and strain) tensors can be derived at any given 
point (x, y); those stress and strain tensors will then automatically satisfy the equilibrium 
equations of linear elasticity Airy stress function is a solution of bi-harmonic equation 

A 2 U = 0. (1) 

It is related to the stress tensor by 

u yy df ' 

o~xx = QyT, (2) 
a - ^ 



X V dxdy ' 

It is well-known j23] that a biharmonic function can be expressed via two analytic 
functions <p(z), x( z ), z = x + iy : 

2U = z<p + zip + x + X- (3) 
It follows from ^ that 

OXX + ^ = 2 + $) , 

o-yy - a xx + 2ia xy = 2 (z& + *) , $ = <f/,V = x" ■ 
Changing variables to polar coordinates (r, 9), we have 

o-„ + o- ee = 2 ($ + $) 
oqb - o„ + 2ia re = 2e 2i0 + *) 

and it follows that 

a„ - ia re = $ + # - e 2ie {W + *) . (6) 
For the displacement u = (u x ,u y ), it follows from the Hook's Law that 

2/i {u x + iu y ) = 2/i (u r + iug) e 10 = —2dU + 2—^ — —(p — K<p — zip — x', 

A + /i (7) 



■2,H , ^ r , J , T , >, (^) 



A + 3/i 



A + ^ 

(equation ([7]) goes back to Lame; the derivation can be found in Ref (22), §32). Lame 
parameters A, /x are related to the Young modulus E and Poisson ratio v as follows, 
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We expand functions $, \I/ in the Laurent series in the outer ring, and in Taylor series 
on the inner disk, 

R^ < \z\ < R? 

(9) 



oo / \ n oo / \ n 



|z| < Ri 

(10) 



oo / \ n / \ n 



n=0 x x n=0 

We assume that the external stress applied to the outer boundary \z\ = R 2 is known, and 
has a Fourier series expansion 

oo 

a rr -ta re = A ne in9 - (11) 

n=— oo 

Thus on the external boundary \z\ = R 2 from ([6]), (|9]), (11) we have that 

(1 -n)a n -Pn-2 + a- n = A n . (12) 

We now turn to analyzing conditions on the internal boundary \z\ = R\. Let us 
introduce dimensionless radius r 1; 

n = §1, (13) 

and let us use the convention that 

b n = q n = 0, n < 0. (14) 

On the internal boundary \z\ = R\, stress components a rr , a r g and the displacement are 
continuous. Using (|6j)-( 10 ) , continuity of such stress components imply that 

(1 - n)a n r\ - p n _ 2 + a„ n r 2 ~ 2n = (1 - n)b n r\ - q n ^ 2 + 6_„r 2 - 2n , (15) 

and the continuity of displacement imply that 

— {{n - l)r 2 a n +p n _ 2 + K 2 a_„rf 2n+2 } = — {(n- l)r\b n + q n _ 2 + kJ)^ 2 ^ 2 } 

A*2 I 1 ! 



Using (12) to eliminate p n , 

_(„_!)(!_ r l) an _ ( n _ i) r 2 &n _ + (i _ ^ + = A „, (17) 

_( n _ 1) (1 _ r l) an + - a _ n (1 + r 2 f 2n K 2 ) - f ± {(n - l)r 2 b n + g n _ 2 + r\- 2n K^ n } 



A, 
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We can solve (17), (18), supplemented with complex conjugates of those equations. We 



denote complex conjugates of coefficients by a, 6, etc. Complex conjugation and n — > —n 



in (17), (18) yields 



(1 



r 2H 


-2n 


' 1 


r ^ 


-2n 


' 1 




A. 


-n- 



) a n + r\ +2n b n - (1 + n) (1 - r?) a_ n + (1 + n)r?6_ n - g_ 2 _ n = A_ n (19) 



From the zero torque condition, we have that 
A = A ; 



(21) 



this is the condition for (17)-(20) with n = to have solution. In fact, solution for n = 



mode is not unique. However, it is clear from ^ that imaginary part of a does not affect 



the stresses, thus we can take 

<2o = do- 

Assuming zero total force, 
Ai = 0; 



(22) 



(23) 



with this condition and using (14), equations (17) , (18) are consistent for n = 1 mode 



and similarly equations (19) , (20) are consistent for n — — 1 



Solving equations (17)- (20), we can find coefficients {ai,bi,qi} ,i e Z : 



a 
bo 



1 2/Ji-/J2+KlA'2 

2 2/ii-/x 2 +Ki^ 2 +'*f (ft2Ail-KlM2+M2-Ml) °' 



1 



Ml+«2Ml 



2 2/xi-^2+KlAt2+»-f (K2/il-Kl£i2+M2-Ml) ' 
Ml+ K 2Ml ^ j 



k _ 

1 fJ-i+K-i/J-2+rf (K2fJ.1-K.1fi2) 

a_i = 0, 



ai 



/Xi+Kl^2 







Qi— n 




b n 




\ <?n-2 / 




?2 > 2, 7"! 



/il+KlM2+»"Y( K 2Ml— K lA*2) 
/ 



n 



l)(l-r?) 



1 



r 



2+2n 
1 



-{n - 1) (1 -? 

\ 1 + r 2+2n K 2 

= Ri 
- R 2 - 



2r™~ 1 + 1 
l + n)(l-r?; 

1 + r\- 2n n 2 
l + n)(l-r?; 



-(n — l)rf 


-1 


\ 


-1 


/ A n 


„2+2n 
'1 









A- n 


-(n-l)rj/i2 


_M2 






A n 


Ml 


Ml 






2 + 2n 

— »1 fti/^2 
Ml 





/ 




\ A_ n 



We note that for the displacement ^ to be single- valued, we must have K!a_ r 

-1 = a_i = 0. 



•a_i 



0, which is satisfied by our solution (24) since a 
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Coefficients {p n } can be found from (12), (24) 



Pn = -M+n ~ (1 + n)a 2+n + O-2-n, 71 G Z. (25) 

We note that it is evident from the above formulas (as well as from dimensional 
analysis) that for a prescribed applied external stress, given as a function of the polar 
angle, the solution for stresses on the domain depends on the ratio of internal and external 
radii r% — jjr, but not the overall scale i?2- 

3. Analyzing solution. 

For a given load, we can compute components of stress tensor using equations ([5]), (|9]), 



(10), (14), (24), (25); those formulas determine stresses created by any transversal load. 

In this paper we will analyze in detail response to uni-directional transversal 
compression applied along a diameter. Such load may be used to estimate a response 
to the transversal compression of a fiber created by impacting projectile, when the fiber 
is pressed against a matrix and other fibers to stay in place. We assume that in such 
case high loads will be applied in a small neighborhood of points of contact, and will be 
uni-directional. Young modulus of the reinforcing outer graphene layer is very high, thus 
we anticipate that material delivering external load (e.g. a projectile and the matrix, or 
loading plate in a static loading experiment) will be deformed at a point of contact with 
the reinforced fiber. To estimate stresses in the Kevlar-graphene system we assume that 
the load is uniform in a small region of contact. 

The load we use in our study is shown in Figure 2. Compressing force is applied on 
the outer boundary in the (-Y) direction in a small neighborhood of Y axes intercept, 
with an opposite balancing force applied in diametrically opposite region. Magnitude of 
the force per transversal to Y surface area element, as a function of the polar angle 8, 

71 

has narrow peaks centered at 9 = ± — ; we take the magnitude to be constant (1 GPa ) at 
such peaks, and we take width of the peak to be w = Outside of the peaks, we take a 
smooth interpolation to a zero value in a small angular region of width e. Thus /y, the 
Y axes projection of the applied force per transversal to Y area element, as a function of 
polar angle is given by 



fy{0) 



i, |0 + f|<^ ( 26 ) 



2 



0, 0±f > 9 w + e. 
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Corresponding stresses applied on the outer boundary , as functions of polar angle, are 
o- rr (0) = f y (0)\ sin0| sin 9, 

e re (6) = ^(6)1 sin 9\ cos 9; 1 ' 

f y as a function of the polar angle 9 is shown in Figure 2 (a). In Figure 2(b) and 2 
(c) we plot components of the stress tensor a rr and a r g, induced in the material by such 
applied load, as a function of the polar angle; here r, 9 are polar coordinates. Stress cr rr 
corresponds to a compression in the vicinity of two diametrically opposite points on the 
Y axes. Stress a r g rapidly changes sign at those points (this is because projection of f y 
onto the polar unit vector eg changes sign). 

To understand how stresses are distributed in the Kevlar fiber and the outer graphene 
layer, we compute components of the stress tensor using equations ([5]), ([9]), (10), (14), 
@, {25}. 

We take experimental values for the transverse Young modulus E\ = 5 GPa and 
Poisson ratio v\ = 0.35 for Kevlar, and E2 = 600 GPa and Poisson ratio i>2 = 0.1 
for the graphene in our analysis, and we use the relationship ^ to estimate the Lame 
parameters A and fi of Kevlar and graphene. For the particular combination of Lame 
parameters which appear in the exact solution (24) we get: /ii = 1.85 GPa, Ki = 1.6 for 
the Kevlar, ^2 = 270 GPa, k>2 = 2.6 for the graphene. We note that using the five elastic 
constants describing anisotropic, transversally isotropic material (when those constants 
are known) will give more accurate determination of the transversal Lame parameters. 
Those five elastic constants were measured for the Kevlar KM2 fiber |25j ; using those, we 
get values of the same order as above for Hi and K\, and the qualitative behavior of the 
solution is similar. 

In Figure 3 we plot the a yy component of the stress tensor on the Y axes (that is on the 
line where the load is applied), for varying values of the ratio of the inner to outer radius 
r\ — : We clearly see that the stress in the Kevlar fiber region is considerably reduced, 
dropping down in the stronger graphene region, with about 2 fold drop when ri ~ 0.96, 
and even bigger drop for smaller r\. The magnitude of the drop varies continuously when 
we vary r\. When r\ approaches 1 , there is no reduction of stress, as expected. 

In Figure 4, we plot all 3 components of the stress tensor, in the Kevlar fiber region, 

for varying thickness of the outer graphene shell. When graphene shell is thin, e.g. 
R 

ri = — ~ 0.99, there is not much reduction in stresses in Kevlar; but for a thicker 
R2 

outer shell, r\ ~ 0.96, there is about 2 fold reduction of stresses in the Kevlar fiber region. 
We see that stresses in the Kevlar fiber region attain a maximum value in the outmost 
region of the fiber and near the Y axes where the load is applied. Stresses in Kevlar fiber 
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region vary continuously with the change of ri, getting smaller when the outer graphene 

shell gets thicker in a continuous manner, as expected. 

In Figure 5, we plot the a rr and ayg components of the stress tensor in the Kevlar 

fiber region, corresponing to polar coordinates r, 9. Those two components of stress tensor 

are selected since those are continuous across the inner boundary between Kevlar and 

graphene regions. Behavior of those components in the Kevlar fiber region is similar to 

the stress components in the Euclidean coordinates (x, y) shown in Figure 4: stresses 

are larger in the outer parts of the fiber and near the Y axes where the external load is 

R 

applied; for a fixed external load stresses decrease continuously with decrease of r 1 = — , 

R-i 

with approximately 2-fold decrease at r\ = 0.96. 



3.1. Estimating yield 

Yield criteria for a broad class of materials may be given in terms of the Von Mises 
stress. The Von Mises stress is proportional to the square root of the sum of squares 
of characteristic values of the trace-free part of the stress tensor. Since Kevlar fiber is 
transversally isotropic, and the transversal yield strength of the Kevlar fiber is much lower 
than the longitudinal one, we will use a 2 dimensional version of the Von Mises stress, 



cm 



yj (o- xx - or yy ) 2 + Aa 2 xy . (2* 



We note that for a planar elasticity problem we consider, the longitudinal stress a z 



g!3 



o, x+/e ) {.°~xx + ffyy) ~ 0.6 [p xx + o-yy), and so is of the same order as planar stresses, while 
longitudinal yield is much higher than the transversal yield. Thus using more accurate 
for anisotropic material yield criteria such as Hill yield criteria will produce similar 
conclusions. 

To illustrate how the Von Mises stress is distributed on the Kevlar-graphene domain, 



we plotted the Von Mises stress (28) on the whole domain, Figure 6 (a), and on the 
Kevlar fiber region, Figure 6 (b), in the case when r% = ^ = 0.95. The external applied 



load is along the Y axes and has a sharp peak of polar angle width ^, as in (26), (27 ). 
Distribution of the von Mises stress is largely similar to the distribution of components 
of the stress tensor in Figure 4. There is a sharp drop of the Von Mises stress occuring in 
the stronger outer graphene layer, with a more than 2 fold decrease of stress in the Kevlar 
fiber compared to the case of no graphene reinforcement. On the Kevlar fiber domain, the 
Von Mises stress is larger in the outmost region and near the Y axes where the external 
stress is applied. 
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Von Mises yield criteria is that yield occurs when Umax , the maximum value of the 
Von Mises stress ou taken over the domain, becomes greater or equal to the critical 
value o-yi e \d , specific for the material. For Kevlar, ciyj e y ~ 2.9 GPa for longitudinal 
loads [10]; but for transversal loads it is lower (e.g. for Kevlar/PEKK composite the 
transversal strength is 21 MPa, |24J). 

If, at a fixed external load, the maximum von Mises stress in the Kevlar region is 
reduced by a factor of 2, say, due to reinforced outer graphene layer, that implies that 
combined Kevlar-graphene system will have twice bigger yield strength (we assume, as 
experimental data suggests, that the yield strength of graphene is much bigger than the 
transversal yield strength of Kevlar, and that linear elasticity approximation suffices for 
the estimate). 

In Figure 7, we show the reduction of the maximum value of the von Mises stress 
in the Kevlar fiber region, as a function of the ratio of the Kevlar fiber radius R\ to the 
total radius R 2 , for various loading cases. 

Figure 7 (a) shows a uni-axial compression in the Y direction, applied to the outer 
graphene layer, and having a sharp peak of angular width ^ on the Y axes. Figure 7(a) 
illustrates typical displacements for such load. In Figure 7(d), we show the reduction of 
the maximum von Mises stress in the Kevlar fiber region, CTmax |^ 2 l as a function of the 
inner to outer radius ratio jj^. 

Figure 7 (b) shows an off-center compression in the Y direction, applied to the outer 
graphene layer at 4 points, having polar angles ±|, Each load has angular width f^. 
Figure 7 (c) shows a one-sided version of such off-center compression, applied at 2 points 
having polar angles ±|. In Figures 7 (b), (c) we illustrate typical displacements for such 
loads. Reduction of the maximum von Mises stress for such loads is shown in Figures 7 
(e) and 7 (f), respectively. 

The calculation reveals that the reduction of the maximum von Mises stress for all 
the three different loading cases was similar. 

Conclusions 

Strengthening of Kevlar-graphene fibers, with the increase of thickness of the outer layer 
of graphene, was predicted by analyzing an explicit solution of the equations of 2D linear 
elasticity. 

Results of computation of von Mises stress in Kevlar-graphene fibers for different 
transversal loadings demonstrated a two-fold increase in the yield strength, when the 
outer graphene shell radius was about 4 % (corresponding to about 8 % weight) of the 
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Kevlar fiber radius. 

The theoretical results presented here can potentially be used to guide experimental 
work, and to model strengthening of other fiber reinforced composites. 
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Figure Captions 

Figure 1. (a) - Schematic diagram of Kevlar-graphene system: the inner part is Kevlar 
fiber with radius R±, the outer layer is graphene with radius R2: (R2 > Ri)', (b) - 
Displacements on the whole domain, in the case g = 0.70 

Figure 2. (a) Applied force per unit transverse area, as a function of polar angle. Force 

71 

is in the Y direction, and is nonzero in an interval of width 9 W — ^ centered at 9 = ± — ; 
(b) o~ rr component of material stress on the outer boundary, as a function of polar angle 
6; (c) a r e component of material stress on the outer boundary, as a function of polar angle 
9 (polar coodinates are denoted r, 9 ) 

Figure 3. a yy component of the stress on the Y axes (that is on the line where the load 
is applied), for varying values of the ratio of the inner to outer radius g : (a) g = 0.90 

(b) g = 0.93 (c) g = 0.96 (d) g = 0.99 Applied force per unit area is the same ,at ±1 
unit, and is in the Y direction; it is applied on the outer boundary at polar angle intervals 
of width centered at 9 — ±|. 

There is a marked drop of stress occurring in the outer reinforcing graphene layer, provided 
that Sa=*i > 4%. 

Figure 4. o- XX) a yy) a xy components of the stress tensor on the Kevlar fiber domain, for 
varying values of the ratio of the inner to outer radius g : 
(a) a xx , § = 0.90 (b) a yy , § = 0.90 (c) a xy , § = 0.90 ' 

(d) a xx , g = 0.93 (e) a yy , f = 0.93 (f) a xy , § = 0.93 
(g) Oxx, f = 0.96 (h) a yy , § = 0.96 (i) a xy , § = 0.96 
(j) **x, f = 0.99 (k) a yy , f = 0.99 (1) a xy , f = 0.99 

Applied on the outer boundary force per unit area is the same, at ±1 unit, and is in the 
Y direction; it is applied on the outer boundary within a polar angle interval of width ^ 
centered at 9 — ±|. 

Figure 5. o rr , a r g polar coordinates components of the stress tensor on the Kevlar fiber 
domain, for varying values of the ratio of the inner to outer radius g : 
(a) <Trv,% = 0.90, (b) <Tre,% = 0.90, 

(c) a„,% = 0.93, (d) are,% = 0.93, 

(e) = 0-96, (f) <Tre,jfc = 0.96, 
(g) a„ , § = 0.99, (h) <Tre,% = 0.99 . 

Applied on the outer boundary force per unit area is the same, at ±1 unit, and is in the 
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Y direction; it is applied on the outer boundary within a polar angle interval of width ^ 
centered at 9 — ±f . 

Figure 6. Von Mises stress o~m (a) on the whole domain , (b) in the fiber region only, in 
the case where ratio of the inner to outer radius ^ is 0.95 . Applied force per unit area 
is ±1 unit, and is a compression along the Y axes; it is applied on the outer boundary 
within a polar angle interval of width ^ centered at 9 = |, with an opposite force applied 
in the diametrically opposite region. 

Figure 7. External loads, corresponding typical displacements, and reductions of the 
maximum Von Mises stress in the Kevlar fiber region as a function of the ratio of the 
Kevlar fiber radius R\ to the total radius R2 : 

(a) uni-axial compression in the Y direction, applied to the outer graphene layer, and 
having a sharp peak of angular width ^ on the Y axes. Typical displacements are shown. 

(b) off-center compression in the Y direction, applied to the outer graphene layer at 4 
points, having polar angles ±|, and of angular width ^ each . 

(c) one-sided off-center compression in the Y direction, applied to the outer graphene 
layer at two points, having polar angles ±|, and of angular width ^ each. 

(d) reduction of the maximum Von Mises stress in the Kevlar fiber region, as a function 
of the ratio of the Kevlar fiber radius R± to the total radius R2, for the central loading 
shown in Fig. 7(a). 

(e) reduction of the maximum Von Mises stress in the Kevlar fiber region, as a function 
of j^, in the case of symmetric off-center loading, Fig. 7(b) 

(f) reduction of the maximum Von Mises stress in the Kevlar fiber region, as a function of 
J^, in the case of one-sided off-center loading, Fig. 7(c). 

Reference red line in (d), (e), (f) corresponds to no graphene reinforcement. 
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Figure 1. (a) - Schematic diagram of Kevlar-graphene system: the inner part is Kevlar fiber 
with radius R%, the outer layer is graphene with radius i?2 : (R2 > Ri)', (b) - Displacements on 
the whole domain, in the case = 0.70 
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Figure 2. (a) Applied force per unit transverse function of polar angle. Force is in 

the Y direction, and is nonzero in an interval of width 6 W = centered at 9 = ± — ; (b) a rr 

component of material stress on the outer boundary, as a function of polar angle 9; (c) a r g com- 
ponent of material stress on the outer boundary, as a function of polar angle 9 (polar coodinates 
are denoted r, 9 ) 
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Figure 3. a yy component of the stress on the Y axes (that is on the line where the load 
is applied) , for varying values of the ratio of the inner to outer radius = 0.90 (b) 

jj 1 = 0.93 (c) ^ = 0.96 (d) jj^ = 0.99 Applied force per unit area is the same ,at ±1 unit, and 
is in the Y direction; it is applied on the outer boundary at polar angle intervals of width ^ 
centered at 9 = 

There is a marked drop of stress occurring in the outer reinforcing graphene layer, provided that 
^>4%. 
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;ure 4. cr xx ,cr yy ,(J xy components of the stress tensor on the Kevlar fiber domain, for varying values 
of the ratio of the inner to outer radius ^ : 
(a) a xx , g = 0.90 (b) a yy , § = 0.90 (c) a xy , g = 0.90 
(d) a xx , § = 0.93 (e) a yy , § = 0.93 (f) a xy , § = 0.93 
(g) ^x, % = 0.96 (h) a yy , f = 0.96 (i) a xy , f = 0.96 
(j) ^x, % = 0.99 (k) a yy , § = 0.99 (1) a xy , § = 0.99 

Applied on the outer boundary force per unit area is the same, at ±1 unit, and is in the Y 
direction; it is applied on the outer boundary within a polar angle interval of width yg centered 
at0 = ±§. 
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Figure 5. a rr , (J r e polar coordinates components of the stress tensor on the Kevlar fiber domain, for 
varying values of the ratio of the inner to outer radius g : 
(a) a rr , § = 0.90, (b) a re , § = 0.90, 
(c) a rr , § = 0.93, (d) a re , g = 0.93, 
(e) a rr , § = 0.96, (f) a r9 , % = 0.96, 
(g) <Trr , § = 0.99, (h) ^ , f = 0.99 . 

Applied on the outer boundary force per unit area is the same, at ±1 unit, and is in the Y 
direction; it is applied on the outer boundary within a polar angle interval of width jg centered 
at<9 = ±§. 
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Figure 6. Von Mises stress ctm (a) on the whole domain , (b) in the fiber region only, in the 
case where ratio of the inner to outer radius jj^ is 0.95 . Applied force per unit area is ±1 
unit, and is a compression along the Y axes; it is applied on the outer boundary within a polar 
angle interval of width ^ centered at 9 = ^, with an opposite force applied in the diametrically 
opposite region. 
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Figure 7. External loads, corresponding typical displacements, and reductions of the maximum 
Von Mises stress in the Kevlar fiber region as a function of the ratio of the Kevlar fiber radius 
R\ to the total radius R2 ■ 

(a) uni-axial compression in the Y direction, applied to the outer graphene layer, and having a 
sharp peak of angular width on the Y axes. Typical displacements are shown. 

(b) off-center compression in the Y direction, applied to the outer graphene layer at 4 points, 
having polar angles ±|, ±^p, and of angular width each . 

(c) one-sided off-center compression in the Y direction, applied to the outer graphene layer at 



Simulation of the Elastic Properties of Reinforced Kevlar-Graphene Composites 



22 



two points, having polar angles ± j, and of angular width yg each. 

(d) reduction of the maximum Von Mises stress in the Kevlar fiber region, as a function of the 
ratio of the Kevlar fiber radius R\ to the total radius i?2, for the central loading shown in Fig. 
7(a). 

(e) reduction of the maximum Von Mises stress in the Kevlar fiber region, as a function of 
in the case of symmetric off-center loading, Fig. 7(b) 

(f) reduction of the maximum Von Mises stress in the Kevlar fiber region, as a function of 
in the case of one-sided off-center loading, Fig. 7(c). 

Reference red line in (d), (e), (f) corresponds to no graphene reinforcement. 



